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Abstract 

Random Dirac fermions in a two-dimensional space are studied numerically. 
We realize them on a square lattice using the 7r-flux model with random 
hopping. The system preserves two symmetries, the time-reversal symmetry 
and the symmetry denoted by {TL, 7} = with a 4 x 4 matrix 7 in an effective 
field theory. Although it belongs to the orthogonal ensemble, the zero-energy 
states do not localize and become critical. The density of states vanishes 
at zero energy as ~ E a and the exponent a changes with strength of the 
randomness, which implies the existence of the critical line. Rapid growth 
of the localization length near zero energy is suggested and the eigenstates 
near zero energy exhibit anomalous behaviour which can be interpreted as 
a critical slowing down in the available finite-size system. The level-spacing 
distributions close to zero energy deviate from both the Wigner surmise and 
the Poissonian, and exhibit critical behaviour which reflects the existence of 
critical states at zero energy. 
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Dirac fermions often appear in condensed matter physics, for example, a transition be- 
tween different quantum Hall states jlfU, two-dimensional graphite sheets a mean- field 
theory of the t-J model ('7r-flux state') and d-wave superconductors 0. It is then nat- 
ural to investigate what happens when disorder is included. Random Dirac fermions in a 
two-dimensional space were investigated by several groups [p|-|T2"|. Possible appearance of 



non- localized states, critical states, in random Dirac fermions was pointed out in 0. Re- 
cently, this disordered critical state was realized in a lattice model, where it was crucial to 
preserve a symmetry denoted by {7i, 7} = with a 4 x 4 matrix 7 in an effective field theory 
@ (see below). 

In this paper, we study the random Dirac fermions numerically beyond the zero modes. 
In order to realize the massless Dirac fermions on a two-dimensional lattice, we use a tight- 
binding model on a square lattice with half a flux quantum (V flux') per plaquette, which 
is described by the Hamiltonian 

Hpure ^} ^ CjtyCj ~\- h.C, (1) 
<i,j> 

where the summation is over the nearest-neighbor bonds. The hopping matrix elements are 
given by tj +Xj j = (— 1)^ and tj +y j = 1, where j = (j x ,jy), % = (1,0) and y = (0, 1). In the 
momentum space, the Hamiltonian is rewritten as 



Hpure - 2 X! ^{kx^y) 

{hx fey ) 



^ cos ky cos k x ^ 



cos k x — cos k y j 
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where the summation is over the magnetic Brillouin zone [— 7r, ix) x [0, 7r) and ipt 



ifex fey ) 



( c \k x k y )i c \k x ky+n))- There are two energy bands E(k) = ±2^Jcos 2 k x + cos 2 k y on the 
magnetic Brillouin zone. They touch at two momenta, k 1 = (k x ,ky) = (tt/2,7t/2) 
and k 2 = (k x ,k 2 ) = (— 7r/2, tt/2). Near the degeneracies k* (i = 1,2), they be- 



have as E(k) « ±2 y (k x — k x ) 2 + (k y — k y ) 2 (i = 1,2). Define continuum variables 
&{x,y) = (ij 1 \x,y),ip 2 \x,y),ip 3 \x,y),ip 4 \x,y)) by cj ~ a[ tf'+hifafa y) +i j ^ip 2 {x, y) + 
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i i x+: > y i/; 3 (x,y) +i ^ ^i/j^x, y) ], where a is the lattice spacing and x = aj x ,y = aj y . Then 
the Hamiltonian becomes in the continuum limit (a — > 0) 
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Therefore our lattice model includes doubled massless Dirac fermions. 

There are several subtleties for the massless Dirac fermions. When the Fermi energy lies 
at zero energy, that is, all the negative energy eigenstates are filled, the Hall conductivity o xy 



is ill defined. An infinitesimal mass determines the a xy in the continuum theory 



Similar 



phenomenon also occurs in a lattice model where an infinitesimal next-nearest-neighbor 
hopping t' opens a gap and the a xy is given by f/\f\ @. Therefore the massless Dirac 
fermions are at a quantum phase-transition point between differnt quantum Hall states. 

Let us consider the effect of randomness in the hopping matrix elements. We set tj +x j = 
(—)■?'» + 6tj +x j and tj +y j = 1 + 6tj +y j, where Stj +X j and Stj +y j are random variables and 
taken at random with constant probability from [-W, W]. It should be noted that this model 
preserves the time-reversal symmetry and belongs to the orthogonal ensemble. Another 
example of the orthogonal ensemble, Dirac fermions with diagonal disorder, was studied in 
PJTD| and it was suggested that all the eigenstates localize, which is consistent with the 



scaling theory of the Anderson localization ||14| . In the case of the random-hopping model, 
however, it was found that the zero mode does not localize but become critical M. Similar 



phenomenon was found at the band center of the quantum Hall states JT^] . In the quantum 
Hall states, however, the time-reversal symmetry is broken and the system belongs to a 
different universality class, the unitary ensemble. In |J, it was confirmed that parameters 
for the critical states form a critical line in the parameter space of the Hamiltonian, which 
is connected to pure massless Dirac fermions at zero energy. We consider that the stability 
of the zero modes against random hopping matrix elements is due to a symmetry of our 
Hamiltonian. The random hopping matrix elements preserves the symmetry in contrast 
to the diagonal disorder. In the language of the lattice Hamiltonian (ID), the symmetry 



means that the transformation c~ 



[—ly^iucj induces sign change of the Hamiltonian. 
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Thus the eigenstates always appear in pairs with energies E and —E. The corresponding 
transformation in the continuum Hamiltonian (^) is given by 7i pm . e — >7^7ip Ure 7 = — 7~t pU re, 
where 7 = <j\ ® <j\ . Since the random hopping matrix elements preserve the symmetry, the 
continuum Hamiltonian 7i for the random-hopping model also satisfies {7i, 7} = 7i r ) + r y'H = 
0. Thus, taking the lowest order in derivatives, we obtain the following form as a possible 
Hamiltonian for the effective field theory 



H = 2ijdx &(x,y) 
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'^{x,y), (4) 



i=l 



where 7 1 = a 2 ® I, 7 2 = <J\ <8> cr 2 , 7 3 = — a 2 <8> en, 7 4 = I <8> cr 2 and Oi(x, y) (i=l, • • -, 4) are 
random variables. 

In this paper, we study random Dirac fermions numerically beyond the zero modes. 
We diagonalize the Hamiltonian for finite squares of size L 2 = 20 2 , 30 2 , 40 2 and 50 2 . To 
obtain reliable statistics, ensemble average over 16000, 16000, 8000 and 3360 realizations 
is performed respectively. The observables are density of states p(E), the dimensionless 
conductance (Thouless number) g(E) and the level-spacing distribution P(s). 

Let us first discuss the density of states p(E) = l/L 2 J2i^(E — Ei). When there is no 
randomness, the p(E) vanishes linearly at zero energy. Recently, whether the density of 
states is finite or not at zero energy for random Dirac fermions is controversial [p|,^2| . The 
p(E) , s for different strength of the randomness are shown in Fig. 1. We have fitted the 
data by the power-law form p(E) = CE a ( w \ Within the numerical accuracy, our results 
support the vanishing density of states at zero energy with an anomalous exponent a(W), 
which depends on strength of the randomness. 

Next, in order to reveal nature of the eigenstates, let us consider the dimensionless 
conductance (Thouless number) g(E). The g(E) is defined by g(E) = V(E)/A(E), where 
V(E) is an energy shift obtained by replacing periodic boundary condition by antiperiodic 
boundary condition and A(E) is a local mean level spacing near the energy E. Numerical 
results for the g(E) are shown in Fig. 2 with L =30, 40 and 50, where ensemble average 
is performed within an energy window whose center is located at each data point. Rapid 



enhancement of the g(E) near zero energy is observed in Fig. 2. It suggests that the 
localization length grows rapidly near zero enrergy. This is consistent with the existence of 
crtitical states at zero energy. One may consider that the zero modes are just on the critical 
point. Then one of the possible scenarios is that the non-zero energy eigenstates are all 
off critical and therefore localized. It suggests an exponential dependence of the finite-size 
dimensionless conductance g(E,L). The g(E,L) obtained numerically decreases when the 
system size increases. It is, however, far from the exponential dependence. In Fig. 3, we 
have plotted the g(E, L) as a function of l/L. It suggests a power-law form g(E, L) oc 1/L 7 
rather than an exponential form g(E,L) oc exp(— L/£). Although we can not exclude a 
possible existence of critical states in a finite energy region, we consider that the non-zero 
energy states may be localized in an infinite-size system and a crossover from the power-law 
form to an exponential form occurs when the system size increases beyond the localization 
length. The localization length of the eigenstates near zero energy may be large compared 
to the available system sizes and we may say that the power-law dependence of the g(E, L) 
is a critical slowing down in the available finite-size system. This also suggests the existence 
of the critical state at zero energy. 

We have also obtained the level-spacing distribution P(s). The P(s)'s of the normalized 
energy separation s = \E n — E n+ i\/ A(E n ) are shown in Fig. 3 and 4, where E n and 
E n+ i are two successive eigenenergies. The P(s) is well described by the Wigner surmise 
P(s) = As^exp(—Bs 2 ) in the metallic regime and becomes the Poissonian P(s) = exp(— s) 
in the insulating regime. The parameter (3 in the Wigner surmise reflects the symmetry 
of the Hamiltonian and /3—1, 2 and 4 for the orthogonal ensemble, the unitary ensemble 
and the symplectic ensemble respectively. The parameters A and B are determined by 
J °° ds P(s) = 1 and / °° ds sP(s) = 1 and, in paticular, A = n/2 and B = 7r/4 for the 
orthogonal ensemble. The P(s) characterizes nature of the eigenstates. States localized in 
different spatial regions are allowed to lie at the same energy. It means that the energy 
levels of the localized states distribute independently, which is described by the Poissonian. 
On the other hand, in metals where the eigenstates are extended, two adjacent energy levels 
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interact strongly, which brings strong energy repulsion and P(s) ~ s 13 near s = 0, where (3 
is determined by the symmetry of the Hamiltonian. The P(s) is well explained by a 2 x 2 
random matrix model, which is the Wigner surmise. Level statistics near the mobilty edge 
has been studied recently and the appearance of critical level statistics is discussed |I6 



Several numerical studies on level statistics near the mobility edge were performed for e.g. 



the three-dimensional Anderson model |18[ and the band center of the quantum Hall states 
19| , which belong to the orthogonal ensemble and the unitary ensemble respectively, and 
it was found that the P(s)'s deviate from both the Wigner surmise and the Poissonian and 
exhibit critical behaviour, i.e., the P(s)'s do follow the Wigner surmise for small s and they 
then deviate from it at higher values of s and show stretched exponential decay. The P(s)'s 
near zero energy for our model are shown in Fig. 3, where the energy window is set [0.1, 0.5] 
| 20fl . The location of the energy window is set sufficiently close to zero energy compared 



to the band width, which corresponds to the energy cut off in the continuum theory. Thus 
we consider that the system is described by random Dirac fermions. We confirmed that the 
P(s)'s near zero energy do not seriously depend on the system size and the location of the 
energy window. We consider that this is due to the fact that, since the location of the energy 
window is set near zero energy compared to strength of the randomness, the localization 
lengths are so long that it exceeds the system size. Thus, although the eigenstates may 
be localized in an infinite-size system, they behave as critical wavefunctions in a finite-size 
system. In fact, the P(s) in Fig. 3 deviate from both the Wigner surmise and the Poissonian, 
and exhibit critical behaviour, i.e., the P(s)'s do follow the Wigner surmise for small s and 
they then deviate from it at higher values of s and show stretched exponential decay. On 
the other hand, when the location of energy window is set in other regions, non-critical 
behaviour is found (see Fig. 4). For example, when L = 20, W = 0.5 and the center of 
the energy window is set E = 1.5, the randomness is so weak that the deviation from pure 
massless Dirac fermion is small and the P(s) is close to the Wigner surmise (see Fig. 4 
(a)), and when L = 50, W = 1.0 and the center of the energy window is set E = 3.75, the 
localization length is within the system size and the Poissonian behaviour is observed (see 
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Fig. 4 (b)). 

In summary, we have studied random Dirac fermions numerically beyond the zero modes. 
Although it belongs to the orthogonal ensemble, the zero-energy states do not localize and 
becomes critical. The density of states p(E) vanishes as ~ E a near zero energy and the 
exponent a depends on strength of the randomness W. It implies that scaling dimensions 
of the operators change with strength of the randomness. It is similar to the case of the 
random gauge-field critical line found in @ . We consider that the existence of the symmetry 
{7i, 7} = is crucial to have the criticality of the zero modes. We have studied nature 
of the eigenstates using the dimensionless conductance (Thouless number) g(E) and the 
level spacing distribution P(s). As is suggested by the numerical results for the g(E), the 
localization length grows near zero energy so rapidly that it exceeds the available system size 
and the observables at non-zero energies exhibits anomalous behaviour, a critical slowing 
down. The P(s)'s near zero energy deviate from both the Wigner surmise and the Poissonian 
and exhibit critical behaviour, as in the case of the quantum Hall states, where rapid growth 
of the localization length and the critical behaviour of the P{s) were observed []I7yi9| . We 
consider that it reflects the existence of critical states at zero energy for our model. As is 
discussed in |7|], the random guage-field critical line is unstable in contrast to the critical line 
of (1+1) -dimensional free bosons, Tomonaga-Luttinger liquid, and a global renormalization- 
group flow for random Dirac fermions is not fully understood yet. Our results show that 
critical states can appear in a lattice model and we consider that our study gives a clue to 
understand a global phase diagram for random Dirac fermions on a lattice. 
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FIGURES 



Fig. 1 The density of states p(E), where L = 50 and W = 0.7, 0.8, 0.9 and 1.0. We have 
fitted the data by the power-law form p(E) = CE a ^ w \ where a(W) = 0.90, 0.74, 0.55 and 
0.39 for W = 0.7, 0.8, 0.9 and 1.0, respectively. We confirmed that the finite-size effect is 
small, comapring the results for L = 50 with those for L = 30 and 40. 

Fig. 2 (a) The dimensionless conductance (Thouless number) g(E), where W = 1.0 and 
L = 30, 40 and 50. (b) l/L-g(E,L) plot for E = 0.30, 0.42, 0.50, 0.62 and 0.70. It suggests 
a power-law form g(E,L) oc 1/L 7 rather than an exponential form g(E,L) oc exp(— L/£) 
in the present system size. We consider that it is a critical slowing down in the available 
finite-size system and a crossover from the power-law form to an exponential form occurs 
when the system size increases beyond the localization length. 

Fig. 3 The level spacing distribution P(s) near zero energy, where L = 50 and ensemble 
average is performed within an energy window [0.1,0.5]. We confirmed that the finite-size 
effect is small, comapring the results for L = 50 with those for L = 20, 30 and 40. We also 
found that there is no substantial difference between results with energy windows [0.1,0.3], 
[0.2,0.4] and [0.3,0.5]. 

Fig. 4 The level spacing P(s), (a) where L = 20, W = 0.5 and ensemble average is 
performed within an energy window [2.0,3.0]; (b) where L = 50, W = 1.0 and ensemble 
average is performed within an energy window [3.5,4.0]. 
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